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Abstract 



The fractional Brownian motion is a generalization of ordinary Brown- 
ian motion, used particularly when long-range dependence is required. Its 
explicit introduction is due to B.B. Mandelbrot and J.W. van Ness (1968) 
as a self-similar Gaussian process W^ H \t) with stationary increments. 

Here self-similarity means that (a~ H W (H) (at) : t > 0) = (W (H \t) : t > 
0), where H £ (0, 1) is the Hurst parameter of fractional Brownian motion. 

F.B. Knight gave a construction of ordinary Brownian motion as a 
limit of simple random walks in 1961. Later his method was simplified by 
P. Revesz (1990) and then by the present author (1996). This approach is 
quite natural and elementary, and as such, can be extended to more gen- 
eral situations. Based on this, here we use moving averages of a suitable 
nested sequence of simple random walks that almost surely uniformly con- 
verge to fractional Brownian motion on compacts when H £ (i, 1). The 

rate of convergence proved in this case is 0(iV- min(ff -4>4)logjV), where 
N is the number of steps used for the approximation. 

If the more accurate (but also more intricate) Komlos, Major, Tusnady 
(1975, 1976) approximation is used instead to embed random walks into 
ordinary Brownian motion, then the same type of moving averages almost 
surely uniformly converge to fractional Brownian motion on compacts for 
any H £ (0, 1). Moreover, the convergence rate is conjectured to be the 
best possible 0{N~ H log AT), though only 0(N~ min ( H '3> log TV) is proved 
here. 
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1 Fractional Brownian motion 



The fractional Brownian motion (fBM) is a generalization of ordinary Brow- 
nian motion (BM) used particularly when long-range dependence is essential. 
Though the history of fBM can be traced back to Kolmogorov (1940)| a nd oth- 



ers, its explicit introduction is due to Mandelbrot and van Ness (1968)| . Their 



intention was to define a self-similar, centered Gaussian process W (H> (t) (t > 0) 
with stationary but not independent increments and with continuous sample 
paths a.s. Here self-similarity means that for any a > 0, 

(a- H W^ (at) : t > o) I (w^ (t) : t > o) , (1) 

where H 6 (0,1) is the Hurst parameter of the fBM and = denotes equality 
in distribution. They showed that these properties characterize fBM. The case 
H = | reduces to ordinary BM with independent increments, while the cases 
H < g and H > \ give negatively, respectively, positively correlated increments, 



see 



Mandelbrot and van Ness (1968)] . It seems that in the applications of fBM 



the case H > ± is the most frequently used. 



Mandelbrot and van Ness (1968)) gave the following explicit representation 



of fBM as a moving average of ordinary, but two-sided BM W(s), s G 
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WiH){t)= Y(H " ^ / d W(s) (t>0), (2) 



where (a;)+ = max(x,0). 

The idea of @ is related to deterministic fractional calculus, which has an 
even longer history than fBM, going back to Liouville, Riemann, and others, see 
Samko et al. (1993)] . Its simplest case is when a continuous function / and 



a positive integer a are given. Then an induction with integration by parts can 
show that 

/«(*) = pTT [\t-sr- 1 f(s)ds 
r (a) Jo 

is the order a iterated antiderivative (or order a integral) of /. On the other 
hand, this integral is well-defined for non-integer positive values of a as well, in 
which case it can be called a fractional integral of /. 
So, heuristically, the main part of (J2J, 
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W a (t) = — / (t-sr- 1 W'(s)ds=— - / (t-sf-'dWis) 
"(a) J Q T{a) J 

is the order a integral of the (in ordinary sense non-existing) white noise process 
W'(t). Thus the fBM W {H \t) can be considered as a stationary-increment 
modification of the fractional integral W a (t) of the white noise process, where 



2 Random walk construction of ordinary Brow- 
nian motion 

It is interesting that a very natural and elementary construction of ordinary 
BM as a limit of random walks (RWs) appeared relatively late. The mathe- 
matical theory of BM began around 1900 with the works of Bachelier, Einstein, 
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Smoluchowski, and others. The first existence construction was given by Wiener 
(1921, 1923) that was followed by several others later. F.B. Knight (1961) in- 
troduced the first construction by random walks that was later simplified by P. 
Revesz (1990). The present author was fortunate enough to hear this version 
of the construction directly from Pal Revesz in a seminar at the Technical Uni- 
versity of Budapest a couple of years before the publication of Revesz's book 
in 1990 and got immediately fascinated by it. The result of an effort to fur- 
ther simplify it appeared in Szabados (1996)] . From now on, the expression 



RW construction will always refer to the version discussed in the latter. It is 
asymptotically equivalent to applying |Skorohod (1965)] embedding to find a 
nested dyadic sequence of RWs in BM, see Theorem 4 in Szabados (1996)] . 



As such, it has some advantages and disadvantages compared to the celebrated 
best possible approximation by BM of partial sums of random variables with 
moment generator function finite around the origin. The latter was obtained 
by Komlos, Major, and Tusnady (1975, 1976), and will be abbreviated KMT 
approximation in the sequel. The main advantages of the RW construction are 
that it is elementary, explicit, uses only past values to construct new ones, easy 
to implement in practice, and very suitable for approximating stochastic inte- 
grals, see Theorem 6 in Szabados (1996)) and also [Szabados (1990)] . Recall 



that the KMT approximation constructs partial sums (e.g. a simple symmetric 
RW) from BM itself (or from an i.i.d. sequence of standard normal random 
variables) by an intricate sequence of conditional quantile transformations. To 
construct any new value it uses to whole sequence (past and future values as 
well). On the other hand, the major weakness of the RW construction is that 
it gives a rate of convergence 0(N~* log AT), while the rate of the KMT ap- 
proximation is the best possible 0(N~i log AT), where is the number of steps 
(terms) considered in the RW. 

In the sequel first the main properties of the above-mentioned RW construc- 
tion are summarized. Then this RW construction is used to define an approxi- 
mation, similar to ([2} , of fBM by moving averages of the RW. The convergence 
and the error of this approximation are discussed next. As a consequence of the 
relatively weaker approximation properties of the RW construction, the conver- 
gence to fBM will be established only for H G (4, 1), and the rate of convergence 
will not be the best possible either. To compensate for this, at the end of the 
paper we discuss the convergence and error properties of a similar construc- 
tion of fBM that uses the KMT approximation instead, which converges for 
all H <E (0, 1) and whose convergence rate can be conjectured to be the best 
possible when approximating fBM by moving averages of RWs. 

The RW construction of BM summarized here is taken from |Szabados (1996)] . 
We start with an infinite matrix of i.i.d. random variables X m (k), 

P {X m (k) = 1} = P {X m (k) = -1} = 1 (to > 0, k > 1), 

defined on the same underlying probability space (O, A,, P). Each row of this 
matrix is a basis of an approximation of BM with a certain dyadic step size 
At = 2~ 2m in time and a corresponding step size Ax = 2~ m in space, illustrated 
by the next table. 

The second step of the construction is twisting. From the independent ran- 
dom walks (i.e. from the rows of Table 1), we want to create dependent ones 
so that after shrinking temporal and spatial step sizes, each consecutive RW 



3 



Table 1: The starting setting for the RW construction of BM 



At 


Aa; 


i.i.d. sequence 


RW 


1 


1 


X (1),X Q (2),X {3),... 


So(n) 


— l^k= 




2~ 2 


2" 1 


X 1 (1),X 1 (2),X 1 {3),... 


Si(n) 


— l^k= 


=1*1 (*) 


2 -4 


2- 2 


X 2 (1),X 2 (2),X 2 {3),... 


S 2 (n) 




=1 X 2 (fc) 



becomes a refinement of the previous one. Since the spatial unit will be halved 
at each consecutive row, we define stopping times by T m (0) = 0, and for k > 0, 

T m (k + 1) = min{?i : n > T m (k), \S m (n) - S m (T m (k))\ = 2} (m > 1). 

These are the random time instants when a RW visits even integers, different 
from the previous one. After shrinking the spatial unit by half, a suitable 
modification of this RW will visit the same integers in the same order as the 
previous RW. (This is what we call a refinement.) We will operate here on each 
point lj G D, of the sample space separately, i.e. we fix a sample path of each 
RW appearing in Table 1. Thus each bridge S m (T m (k + 1)) — S m (T m (k)) has to 
mimic the corresponding step X m _i (k+1) of the previous RW. We define twisted 
RWs S m recursively for m = 1, 2, 3, . . . using S m —i, starting with So(n) = So(n) 
(n > 0). With each fixed m we proceed for k — 0, 1, 2, . . . successively, and for 
every n in the corresponding bridge, T m (k) < n < T m (k +1). Any bridge is 
flipped if its sign differs from the desired (Figs. 1-3): 

X , n) = { X m (n) if S m (T m (k + 1)) - S m (T m (k)) = 2X m _ 1 (k + 1), 
— X m (n) otherwise, 

and then S m (n) = S m {n—1)+X m (n). Then each S m {n) (n > 0) is still a simple, 
symmetric RW, see Lemma 1 in Szabado s (1996)] . Moreover, the twisted RWs 
have the desired refinement property: 

^S m {T m (k)) = S m -i(k) (m > 1, k > 0). (3) 




Figure 1: B Q (t;u) = S {t;uj). 

The last step of the RW construction is shrinking. The sample paths of 
S m (n) (n > 0) can be extended to continuous functions by linear interpolation, 
this way one gets S m (t) (t > 0) for real t. Then we define the mth approximation 
of BM (see Fig. 4) by 

B m (t) - 2- m S m (t2 2m ). (4) 
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X 




Figure 4: Bi(t;u>). 



Compare three steps of a sample path of the first approximation Bo(t; uj) and 
the corresponding part of the second approximation Bi(t;u) on Figs. 1 and 4. 
The second visits the same integers (different from the previous one) in the same 
order as the first, so mimics the first, but the corresponding time instants differ 
in general: 2~ 2 Ti(fc) ^ k. Similarly, ([3]) implies the general refinement property 

B m+1 (T m +i(k) 2- 2 (" l+1 >) - B m (k2- 2m ) (to > 0, k > 0), (5) 
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but there is a time lag 

T m+1 (k) 2- 2{m+1 '> - k2- 2m ^ (6) 

in general. 

The basic idea of the RW construction of BM is that these time lags become 
uniformly small if m gets large enough. It can be proved by the following simple 
lemma. 

Lemma 1. Suppose that X\,X2 1 ■ ■ ■ , Xjv is an i.i.d. sequence of random vari- 
ables, E(Xk) — 0, Var(JTfc) = 1. and their moment generating function E (e uXk ^ 
is finite for \u\ < uq, uq > 0. Let Sj = X± + • • • + Xj, 1 < j < N . Then for 
any C > 1 and N > Nq(C) one has 

max ISA > (2CN\ogN)^ \ < 2N 1 ~ C '. 

0<j<N J J 

This basic fact follows from a large deviation inequality, see e.g. Section 
XVI, 6 in Feller (1966)] . Lemma [T] easily implies the uniform smallness of time 
lags in 

Lemma 2. For any K > 0, C > 1, and for any m > mo(C), we have 

P | o<fc max < ^ |T m+ i(fc) 2- 2 ( m+1 ) - k2~ 2m \ > Qc^Clog, k) ' 2-™| 

< 2(K2 2m ) 1 - c , 

where log„(a;) = max(l, log x). 

Not surprisingly, this and the refinement property ([5]) imply the uniform 
closeness of two consecutive approximations of BM if m is large enough. 

Lemma 3. For any K > 0, C > 1, and m > m\(C), we have 

max \B m+1 (k2- 2m ) - B m (k2- 2m )\ > i^flog* K)im2-^ 

0<k2~^<K 

< iim 2 ™) 1 - . 

This lemma ensures the a.s. uniform convergence of the RW approximations 
on compact intervals and it is clear that the limit process is the Wiener process 
(BM) with continuous sample paths almost surely. 

Theorem 1. The RW approximation B m (t) (t > 0, m = 0, 1, 2, . . .) a.s. uni- 
formly converges to a Wiener process W{t) {t > 0) on any compact interval 
[0, K],K > 0. For any K > 0, C > 3/2, and for any m > 7712(C); we have 

max K \W(t) - B m (t)\ > K i(log* K)im2-^ < §(K2 2m ) x - c . 

The results quoted above correspond to Lemmas 2-4 and Theorem 3 in 
Szabados (1996)] . We mention that the statements presented here are given 



in somewhat sharper forms, but they can be read easily from the proofs in the 
above reference. 



G 



3 A pathwise approximation of fractional Brow- 
nian motion 



An almost surely convergent pathwise construction of fBM was given by Car- 
mona and Coutin (1998) representing fBM as a linear functional of an infinite 
dimensional Gaussian process. Another pathwise construction was given by De- 
creusefond and Ustiinel (1998, 1999) which converges in the L 2 sense. This 
construction uses discrete approximations of the moving average representation 
of fBM @, based on deterministic partitions of the time axis. More exactly, ^ 
is substituted by an integral over the compact interval [0,i], but with a more 
complicated kernel containing a hypergeometric function too. 

The approximation of fBM discussed here will also be a discrete version of 
the moving average representation @ of fBM, but dyadic partitions are taken 
on the spatial axis of BM and so one gets random partitions on the time axis. 
This is asymptotically a Skorohod-type embedding of nested RWs into BM. As 
a result, instead of integral we have sum, and BM is substituted by the nested, 
refining sequence of its RW approximations discussed in the previous section. 
Since ^ contains two-sided BM, we need two such sequences: one for the right 
and one for the left half-axis. From now on, we are going to use the following 
notations: m > is an integer, At = 2~ 2m , t x — xAt (x g R). Then by 



definition, the rath approximation of fBM is: B, 
integers k, 



(H) 



(0) = 0, and for positive 



1 fe_1 i 

• XV E [( tk ~ _ Hr)+~*] [ B m{U + At) - B m {t r ) 



T{H+\) 



r— — OQ 



-2Hm 



k-1 



T (H 



n E [(k-r) H -i-(-r)^]x m (r + l), (7) 



r— — oo 



where the convention H 2 = is applied even for negative exponents. 
B^n^ is well-defined, since the "infinite part" 



—1 — 1 

E [(k-r) H -i-(-r) H -i]x m (r+l)=: £ Y v 



'k,-r 



converges a.s. to a random variable Zk by Kolmogorov's "three-series theorem" : 
E(lft,„) = and 



k 

1 + - 

v 



H- 



-, 2 



- 1 



E 



const 
7,3 — in 



< 00. 



( H) 

It is useful to write B m in another form applying a discrete version of 
integration by parts. Starting with ([7]) and rearranging it according to B m (t r ), 
one obtains for k > 1 that 



B%\tu)= E 



h(t r - At,t k ) - h(t r ,t k ) 
At 



B m (t r )At, 



(8) 
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where we introduced the kernel 
h(s,t)= 1 



r(H + i) 

This way we have got a discrete version of 



(s < t). 



(9) 



W {H) {t) = 



-1 



r(H + i) 



{t-s) H -? -(-«)+ 5 W(s)ds, (10) 



which is what one obtains from ([2]) using a formal integration by parts (cf. 
Lemma [5] below) . 

To support the above definition we show that B^ has properties analogous 
to the characterizing properties of fBM in a discrete setting. 

(a) Bm is centered (clear from its definition) and has stationary increments: 
if ko and k are non- negative integers, then (substituting u = r — k^) 

B^\t ko +t k )-B^\t ko ) 

0-2Hm (ko+k-l 



r h 



TVi 12 (ko + k-r) H -?X m (r + l) 



r=0 



J2(ko-r) H -?X m (r + l) 
r=0 

53 [(ko + k-r) H -? -(k -r) H -?]x m (r + l) 



r— — oo 



-2Hm 



r h 



T\\ 12 (k-u) H -2X m (k + u + l) 



>u=— fen 



Hi? 



53 (^-^(fco + n + l) 

=— fe 

53 [(* - u) ff ~* - (~u) H -i\ X m (k +u + l)\ 

jy ^ [(fc-u)^ -(-«)+-*] X m (*0 + « + l) = B^(t*). 



li — — oo 



-2Hm 



U— — OQ 



(b) Bffi is approximately self-similar in the following sense. If a = 2 2m °, 
where mo is an integer, mo > — m, then for any k non-negative integer for which 
ka is also an integer one has that 



a- H B^\ak2- 2m ) 

■p- ^ ^ak-r) H -i-(-r)^]x m (r + l) 



a~ H 2- 



r— — oc 



2~ 2i? (m+mo) 



i 5^ mo (fc2- 2 ™). 



TT" E [(fc2 2roo -r) ff -i-(-r)^-*]x ro (r + l) 
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On the other hand, Lemma|3](and Theorem^]) below show that Bm and B^K 

( M ) 

(and B m ^ n ) are uniformly close with arbitrary large probability on any compact 
interval if m is large enough (when H > i). It could be proved in a similar 
fashion that for a = j, where j > is an arbitrary integer, 2 2n < j < 2 2 (™ +1 ) 
with an integer n > 0, the finite dimensional distributions of 

n-H(2m+log 2 j) ok-i 

a- H BW(ak2-* m )= 2 v £ [<Jk - r)*~* - (-r)J~*] X m (r+1) 

V 2/ r— — oo 

can be made arbitrarily close to the the finite dimensional distributions of B m J_ n 

( M\ 

if m is large enough. Consequently, Bm is arbitrarily close to self-similar for 
any dyadic a = j2 2m ° if m is large enough. 

(c) For any < t± < ■ ■ ■ < t n , the limit distribution of the vector 

(B^\t^B^\tt ] ),...,B^(t^)) 

as m — > oo is Gaussian, where t { ™ ] = [t.j2 2m \2~ 2m , 1 < j < n. This fact follows 

from Theorem [2] (based on Lemma [5]) below that states that the process Bm 
almost surely converges to the Gaussian process on compact intervals. 



4 Convergence of the approximation to fBM 

At first it will be shown that two consecutive approximations of fBM defined by 
([7]), or equivalently by ([5]), are uniformly close if m is large enough, supposing 
H > j. Apparently, the above RW approximation of BM is not good enough 
to have convergence for H < j. 

When proving convergence, a large deviation inequality similar to Lemma [1] 
will play an important role. If X\, X%, ... is a sequence of i.i.d. random vari- 
ables, P {Xk — ±1} = h, and S — J2 r o, r X r , where not all s r £l are zero and 
Var(S) = J2 r a l < °°> then 

P{\S\ > x(Var(S*))^} < 2e"T (x > 0), (11) 

see e.g. p. 33 in |Stroock (19 93) . The summation above may extend either to 
finitely many or to countably many terms. 

As a corollary, if Si, S2, . . . , Sn are arbitrary sums of the above type, one 
can get the following analog of Lemma [TJ For any C > 1 and N > 1, 

pj max \S k \> (2C log N)i max (Var(S fe ))H 

[l<fc<AT l<fc<W J 

N 

< ^p{|S* fc | > (2Clog^Var(S* fc ))5| < 2Ne- cl ° sN = 2N 1 - . (12) 
k=l 

Lemma 4. For any H £ (i 1), K > 0, C > 3, and m > 1713(C), we have 
P | o< max x iB^tk) - B%\t k )\ > a(H, K)m2^ m J < %{K2 2m ) l - c 
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where tk — k2 2m for k > integers, f3(H) = min(2if — i, ~) and 



a(H,K) 



T(H+l) 



(l-ff)l 



\8K* + 36 


ff-i 


AT*"*) 




2 





a(H,K) = 



(log, jQj 

r(F + i) 



Iff — -I 

-^ + (log,fr)t (5 + 312 



(1-ff)* 



ff- - 
2 



fv 



ff-i 



if H £ l). fTTie case ff = | is described by Lemma\^) 

Proof. The proof is long, but elementary. Introduce the following abbreviations: 
AB m (t) — B m (t + At) — B m (t), AB m+ i(t) = B m+ i {t+\At)-B m+1 {t). Using 
(J?]) and then substituting u — Ar + j, one gets that 

BS 1 (t fe ) = sS 1 (4fc2- 2 (-+ 1 )) 

ry E [(4fc-u) H -*-(-u)*- 3 ]x m+1 (u + i) 



t(h + 


1) 




-1 


r(ff + 


1) 




(4r 


1 




rfff + i) 



U— — 00 
fc-1 3 

E E 

r— — 00 j— 

+ i + i) 

fe-i 3 

E E 

r— — 00 j=0 



k-r- J - 
4 



H- 



AB m+1 {t r 



So, subtracting and adding a suitable "intermediate" term, one arrives at 



fc-1 3 



E E [(**-*h-4) h_ * - Hh-*)? - *] As m+ i(w 4 ) 



r— — 00 j=0 



fc-1 3 



ff- 



E E{[^-W^-(-W+ 

r— — 00 j=0 

- [(i* - t r ) H ~i - (-tr)J - *] } AB m+1 {t r+i ) 

it n r 



fc-1 3 



7' — — OO j—0 

{Z m ,k + Ym,k + V rn .k + f/ m ,ifc)- 



-Af? m (i r ) 



(13) 
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Here we introduced the following notations: 



fc-1 3 

Zm.k = EE [^ k _ t r+i) H ' i - (** " *r) 
fe-l 3 



r=0 j=0 



(fc-r- _ (fc-r) ff -3 



A_B m+ i(£ r+ i) 

X m+1 (4r+i + l) (14) 



and 



fc-i 



3 r 



r=0 
fc-1 



j=0 



1 



AB m+ i(i r+ i) — -AB m (t r ) 



Efa — *r)' H 2 {[-B m+ i(t r+ i) — -B m+ i(i,-)] — [B m (t r+ i) — B m (t r )]} 



r=0 
k 



= ^[{tk-t r -i) H -^ -{t h -t r ) H -^ [B m+1 {t r ) - B m (t r )] , (15) 
r=0 

applying "summation by parts" in the last row, as in Similarly, we in- 
troduced the following notations for the corresponding "infinite parts" in (fT5|) 
(using v = — r): 



V, 



= 2 
and 

Um.k 



-2Hm- 



oo 3 

'EE 

v=l j=0 



(k + v- J -) H -^ - (v - J -) H - L * — (k + v) H - L * + v H -i 



xX m+1 (-Av + j + l), (16) 



oo o 

Eh+^-** _i ]E 



v=l 



AB m+ i(— t v _i) — -AB m (—t v ) 



oo 



v=l 



{[B m +i(— t v -i) — B m+ i(— t v )] — [B m (— t v -i) — B m (—t v )]} 



oo 

= £[(' 



fe + t v+1 ) H ~i - t"} - (t fc + t v ) H ~i + t" { 



v=l 



x{B m+1 (-t v )-B m (-t v )]. (17) 
The maxima of Z mt k,Y m ^,V mt k and U m .k can be estimated separately: 



(18) 



^ T777T — iT ( max l-^m.fcl + max l^m.fcl + max |V m ,fc| + max | C7, 

1 [H + ^ j \ k k k k 



m,fc| I : 



where each maximum on the right hand side is taken for 1 < k < K2 2m and 
one can suppose that K2 2m > 1, that is, At < K, since otherwise the maximal 
difference in (fH?)) is zero. 
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(a) The maximum of Z m ^. In the present case the large deviation inequality 
(fTTT) . or rather, its corollary (|T2"j) is applied. By (fT¥|) . 



k—l 3 



Var(Z m , fe ) = 2 



-AHm-2 ' 



EE 

r=0 J=0 



~i 2 



(1 



4(fc-r)' 



3 - 1 



The term in brackets can be estimated using a binomial series with < j < 3, 
k-r > 1: 



4(fc- 



< 



Thus 



ii 



°° /TT 1 

E( 5 )< " 

-1 



4(fc - r) \ 4(fc - r) 



< 



1 



4(fc - r) 
J 



4(fc-r) 



1 



3 

E 

J=0 



1 



4(/c - r) 
Also, if fc > 1, one has 



< H 



1\ 91 1 



9 {k-r) 2 



E( fc 

r=0 



\2ff-3 



fc-1 



< 1+ / (fc-^^-^dx < 



3 1 



2 1-F 



Then for any k > 1 it follows that 



Var(Z m , fc ) < 2 



-4flm 



i? 



1\ 273 1 



2 / 72 1 - i? 



Hence taking TV = if2 2m and C > 1 in (fl2"|). one obtains that 

(1 - H)-h- 2Hm (2ClogN)^ 

< £p{|Z m , fe | > (2ClogiVVar(Z m , fc ))*} < 2jV 1 " c . 



l</c<JV 
JV 



/273\ 2 








V 72"J 





fc=l 



Since 

logiV = log (K2 2m ) < (1 + log4)mlog„ if < 2.5m log, if, 
one obtains the following result: 



max \Z m k \ < 5 

Kk<K2 2m 



H - - 

2 



(1 - H)~*(C log* K) i to32 



-2Hm 



(19) 



(20) 



with the exception of a set of probability at most 2 (if 2 
if > and C > 1 are arbitrary. 



-"A 1 c w here m > 1, 
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(b) The maximum ofY m> k- By its definition (|15l) . 



max |y m ,fe| < max \B m+1 (t r ) - B m (t r )\ 

Kk<K2 2 ™ 0<f r <Jf 



x max V \{t k -t r -i) H i-(t k -t r ) H 2 
a<t k <K ^— ' l v 
o<t r <t fc 



The first factor, the maximal difference between two consecutive approximations 
of BM appearing here can be estimated by Lemma [3] For the second factor one 
can apply a binomial series: 



r=0 



|(*fc - tr-l)*"* - (tfc - tr)*~* 
2 -m(2H-l) | x 

2— m(2H-l) ) ^ 



< 2 -m(2H-l) J 1 



yH-i 



fe-2 



r=0 
fe-2 

i| + E( fc 

r=0 
fc-2 

l \ + Y t {k-r) H 



k- 



E 



H 



1 \H — 7j\ 
' 2 J ±1 

k — r 



{k - r) s 



Since for ii ^ | 

fc-2 
r=0 



2 < 



fc-1 



(fc-l) H ~5dl 



l-k 



H- 



it follows for any m > that 

fc 



max 

Kfc<K2 2 



< 2- m < 2H - 1 ) max (l 

Kk<K2 2m I 



1 



1 - k 



H- 







-2m{H- 


i) 




' 1 






i) 


K H ~i l) 




I 2 - 


-2m(H- 



if < i? < |, 



(In the last row we used that here 2 2m < JsT.) 

Combining this with Lemma [3l we obtain the result 



max \Y m k \ < 

Kk<K2 2m 



3Ki(\og,K)im2- 2m ( H -^ if \ < H < ±, 
2 H -^K H -i{\og^K)Tm2-^ if § < H < 1. 



(21) 



with the exception of a set of probability at most 3 (K2 2m ) 1 C , where if > 0, 
C > 1 are arbitrary, and m > m\{C). Thus in the case < H < | we have only 
a partial result: the relative weakness of the above-described RW approximation 
of BM causes that apparently we have no convergence for < H < i 
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(c) The maximum of V m ,k- Here one can use the same idea as in part (a), 
including the application of the corollary ([12]) of the large deviation principle. 
We begin with (TTtjI) . 

Var(F m , fc ) 



oo 3 

_ 2~ 4fjm ~ 2 

v=l j=0 
oo 3 ( 

v=l j=0 I 



(k + v — 3 -) H -i -(k + v) H - L * -(v- ] -) H -~i + v H -^ 



4(fc + v) 



H-i 



- v H ~h 



1 



J_ 

4c 



H-i 



1 



As in (a), now we use binomial series for the expressions in brackets (k > 1, 
< j < 3, w > 1): 



A = (k + v 



H-i 



1 - 



J 



4(fc + v) 



H-i 



s=l 



H 



- 1 
j 



4(fc + v) 



and 



1-f 

Av 



H-i 



Then A and _B have the same sign and < Ai < |^4| < |£?| < B%, where 



Ai = (fc + «) 



4(fc + v) 



H ~2 



and 



Hence 



B 2 = v 



1 

2 






K) 


-l 

< 






*-5 


4w 






'"5 



4-j 



(A — B) 2 < A\ + B 2 < [H 



Since for any A; > 0, 



(fe + v) 



2H-3 



4-j 



it follows that 



^(fc + «) 2ff - 3 < 1 + / (k + x) 2H - 3 dx = l + 



(k + 1) 



2H-2 



< 



,2H-3 



3 1 



2-2H 21 — H 



Var(V m , fe ) 



< 2 



-AHm-2 



H 



2 oo 3 

EE 

v=lj=0 



{k + vf H - A + 



4-j 



„2H-3 



< 



791 (#-2-) 2 _ 4gm 
192 1 - if 
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Applying corollary ([12"]) of the large deviation inequality with TV = K2 2m one 
obtains that 



max |F m ,fc| 

Kk<N 



> 



/79l\ 2 








V192J 





< 



(1 - H)-*2- 2Hm (2C\ogN)2 
J2 P {\ V ™,k\ > (2ClogArVar(F m , fe ))'} < 2N 



l-C 



k=l 



Hence using (fT9")l one gets the result 



max \V m ,k\ < 5 

Kk<K2 2m 



H 



(1- if) " 2 (C log, iT) 5 m i2 



-2Hm 



(22) 



with the exception of a set of probability at most 2 (K2 2m ^ 1 C , where m > 1, 
if > and C > 1 are arbitrary. 

('dj TTie maximum of U m ,k- We divide the half line into intervals of length 
L, where L > 4K. For definiteness, choose L — AK . Apart from this, this part 
will be similar to part (b). In the sequel we use the convention that when the 
lower limit of a summation is a real number x, the summation starts at \x], 
and similarly, if the upper limit is y, the summation ends at [y\ . By (fT7|) . 



\U, 



in, k 



^E E 

3=1 (j-l)L<t v <jL 



,H- 



(tk + t«+l) 3 - (t k + t v )"-2 - t v + 1 2 + t 

x \B m+ i(—t v ) — B m (— t v )\ 

00 

< V max \B m+1 (-t v ) - B m (-t v )\ (At) H "3 

' U-l)L<U<iL 



3 = 1 



(k + v + - {k + v) H ~^ - (v + 1) H ^ + v H ~? 

u=(j-l)L2 2m + l 



(23) 



Lemma [3] gives an upper bound for the maximal difference between two 
consecutive approximations of BM if j > 1 is an arbitrary fixed value: 



max \B m+1 (-*„)- B m (-t v ) 

U-l)L<t v <jL 

< (jL)i (log,(j^)) l m2-T 



< 



L« (log,, L) 4m2 ™ if j = 1 , 

2j*(log* j)ii3(l g, L)37772"^ if j > 2, 



(24) 



(25) 



with the exception of a set of probability at most 3 (jL2 2m ) 1 ° , where C > 1 
is arbitrary and m > mi(C). This implies for any C > 3 and m > TOi(C) that 
the above inequality (|2"5)l holds simultaneously for all j = 1, 2, 3, . . . with the 
exception of a set of probability at most 



3(L2 2m )^ C J2j^ C < 3(L2 2m )^ C — < (K2 2m f- C . 



(26) 
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For the other major factor in (|23p binomial series are applied as above, with 
m > 0, k > 1, and v > 1: 



A = (fc + v + 1) H ~* - (fc + v) w ~* = (fc + u) 



i 



1 



k + v 



H- 



v ; ^\ s J(k + v) s 

s—1 v ' v ' 



and for v > 2: 



b = o+ir~5_ w 



H- 

1 + -I -I 



= « H -*E 



h - i\ i 



while B — 2 H 2 — 1 when w = 1. Then A and -B have the same sign, < A\ < 
\A\ < \B\ < B 2 , and so \ A-B\<B 2 -A 1 , where 



Ai =(k + v) 



-(k + v) 



1 










(H 


2 







andB 2 =w H -i|i/-i| . 

Thus if the second major factor in (|23p is denoted by C m ^k,j, we obtain for 
any j > 1 that 



3 L2 2 



(Ai) ff -2 y (k + v + l) H -i -(k + v) H ~i -(v + l) H ~i +V 1 

v={j-l)L2 2m + l 



E {v H -i-(k + v) H -i + \ (l- H ){k + v) B -i\ 



jL2 



D=(j-l)L2 2m + l 

For H ^ h one can get the estimates for j = 1: 



L2- 



E < 1 



L2" 



c ff 2 dx = 



(At) 



i i? 



l ' 



and for j > 2: 

jX2" 
I 

u=(j-l)L2 2 ™+l 



W"" ,jL2 2 



7T_ ^ 

x 2 ax 



0-l)L2 2 " 



(At) 



r-H 



UL) H -*-((j-l)L)"-* 



H- 
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further, for any j > 1, 



3^ Am ,jL2 2m + l 

E {k + v) H -i> {k + xf-idx 

v=(j-l)L2 2 ™+l 



'(j-l)L2^+l 



and also for any j > 1, 

jL2 2 



(tfc+i + _ ( tfc+1 + (j _ i)L) ff 3 



E (k + v) H -i< 



«=0-l)L2 2m + l 



0'-l)£2 2 



-if 



(t fc + (j-l)i) H -^-(t fe +. 7 L) ff -I 



Denote the sign of a real number x by e x (0 if x = 0). When j = 1, it follows 
that 



Cm,fc,i < e H-i { L H 



1-1 



ifc+i 



1 + 



if- 



+ ( # - 2 ) ( A *) 



+ 



H 



At 



L H-i 



i + 



tk — L 



1 + 



tk 



H- 



and similarly, when j > 2, 

C m , fe ,, < e ff i [(ji)* - * - ((j - 1)L) H ^ 



-(tfc+i +jL) H -? + {tk+i + {.] - l)L) 



H-l 



+ 





At r 


H- 1 - 


2 


2 I 



(t k + {j-l)Lf-i -( tk +jL) H -i 



= e H _i{jL) 



H-l 



1- 1-- 



H- 



1 + 



tk+1 



H- 



+ 1 + 



+ 



At 



1 + 



tk — L 



H-i 



1 + 



Applying binomial series here again, first we get when j > 2 that 



-H- 



i-li-i 



if- 



e ^E 



< 



< 



ff - - 

2 



H ~2 



1 1 

7 + 5 



if - - 

2 



"4 
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since each term of the series is positive. Furthermore, with any j > 1, 



-H- 



L — tk+i 



H- 



E 



JL 

H-h\ (-1)' 



JL 



H- 



[(L-t k+1 y-(-t k+1 y] < 



H ~~2 



since each term of the series is negative: L = 4K > 2t k +i, and the term in 
brackets is not larger than 2(L — At) s . Finally, 



1 - 



L_tk 



H- 



if- 



E 



(-1) 6 



[(L-t k y-(-t k y 



L- At 
JL 



4L 



since each term of the series is positive and the term in brackets is not larger 
than §(Z - At) s . 

Thus when j > 2 it follows for any m > 0, fc > 1 that 



H — 



- H 



r 



H- 



At 



4L 



< Off 



2 



where a#- 



if < H < |, 
if | < if < 1. 



In a similar manner, when j = 1 one can get for any m > 0, fc > 1 that 



C m ,i;,l < £ff_li i 



2 


1 

2 — 




1 

~ 2 




2 


—L H 
2 


■'( 


H 



-c H _i^- 5 ](At)^ 

AL 
3At 



2 

3 HL 



(At) 



H-i 



l(At) H -^ if < if < i. 



if | < H < 1. 

Then combine these results with (|2"5j) and (j2"6")l in (j2"5)l . Using 

OO „QQ 

El , .3 ri_5 / rj-_ 9 

J 4 (log*j) 4 j 2 < / a; ilog^da; 



J=2 



/e /' oc 

x 3 dx + / a; 3 log cc da; 

-l 



- H 



-H 



H- 



< 



2.5 if0<#<i, 
16.5 ifi<i?<l, 



(27) 
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one can get the result of part (d). Consider first the case \ < H < 1: 



max \U mk \ 

Kk<K2 2m 



< LiQog,L)-*m2-r-L 



H- 



< 









^3 + 312 




) 









'if— 



33 L*(log^L)3 TO 2"T 



4 (log,, K) l m2~ 



H - 



-H- 



(28) 



for any C > 3 and m > mi(C) with the exception of a set of probability at 
most (K2 2m ) 1 ~ C . (Recall that L = AK.) 

In the second case when < H < \ the above method apparently gives 
convergence here (just like in part (b)) only when \ < H < |: 

max |C/ m ,fe| 

Kk<K2 2m 



< Li(log,L)im2-^-(At) H -i +5 L*(log„ L)*m2~T - 



< 5if3(log^)lm2- 2m ( H -3) +36 



^1^,^)^2-7, (29) 



for any C > 3 and to > mi(C) with the exception of a set of probability at 
most (K2 2m ) 1 ^ C . 

Now one can combine the results of parts (a), (b), (c), and (d), see (fH?]l . 
([20]) . (|2TJ) , p2]l. (|28|. ([29)1 . to obtain the statement of the lemma. Remember 
that the rate of convergence in parts (a) and (c) is faster than the one in parts 
(b) and (d) . Particularly, observe that there is a factor to in (b) and (d) which 
has a counterpart m5 in (a) and (c). Since in the statement of this lemma we 
simply replaced the faster converging factors by the slower converging ones, the 
constant multipliers in (a) and (c) can be ignored if to is large enough. □ □ 

It is simple to extend formula ((8]) of the rath approximation Bm of fBM 
to real arguments t by linear interpolation, just like in the case of the rath 
approximation B m (t) of ordinary BM, see e.g. in Szabados (1996)] . So let 
to > and k > be integers, 7 € [0, 1], and define that 

B^(tk +J ) = 7i4 H W) + (1 - i)Bg Htk) 
k 



(30) 



+ [(-*r) 



TIT £ [(** 



(tk £ r ) 2 B rn it r j r ^ 



2/ r— — 00 
H 



- (-*r-l) 



Bmi^r) • 



Then the resulting continuous parameter approximations of fBM Bm (t) (t > 0) 
have continuous, piecewise linear sample paths. With this definition we are 
ready to state a main result of this paper. 

Theorem 2. For any H e (±, 1), the sequence B^\t) (t > 0, m = 0, 1, 2, . . . ) 
a.s. uniformly converges to a fBM W^(t) (t > 0) on any compact interval 
[0,K], K > 0. If K > 0, C >3, andm> m 4 (C), it follows that 



0<t<K 



\wW(t)-B£\t)\> 



a(H, K) 



r TO2-' 3 ^ m I < g^ 2 " 1 ) 1 - , 
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where a(H, K) and (3(H) = min(2i? — ^, i) are the same as in Lemma^ (The 
case H = k is described by Theorem^) 

Proof. At first we consider the maximum of IB^/^t) — Bm '(t)\ for real t £ 
[0, if]. Lemma S] gives an upper bound D m for their maximal difference at 
vertices with t = t k = kAt: 

max |B<3i(**) - B *n >(**)! < D m , 

except for an event of probability at most 8 (if 2 2m ) 1 ~ c . Since both i?£f+i(i) 

and Bin (t) have piecewise linear sample paths, their maximal difference must 
occur at vertices of the sample paths. Let M m denote the maximal increase of 
Bi n between pairs of points tk,tk+i in [0, if]: 

max \B^(t k+1 ) - B^\t k )\ < M m , 

U S. * k S " 

except for an event of probability at most 2(K2 2rn ) 1 ~ c \ cf. (f!I2j) below. A 

sample path of 5^+1 W ma kes four steps on any interval [tk, tk+i]- To compute 
its maximal deviation from D m it is enough to estimate its change between the 
midpoint and an endpoint of such an interval, at two steps from both the left 
and right cndpoints: 

max {B^t^i) - B^(t k )\ < 2M m+1 , 



except for an event of probability at most 2(K2 2 ^ n+1 ^) 1 c . Hence 
max, \B^ (t k+i ) - S W (t k+i ) \ 



0<t k <K 



= max 

0<t k <K 

< max \B£l(t k ) - B^(t k )\ + max \B^(t k±i ) B^t^ 

< D m + 2M m+ i, 

except for an event of probability at most (8 + 2 3 ~ 2C )(K2 2m ) 1 ~ c . The expla- 
nation above shows that at the same time this gives the upper bound we were 
looking for: 

max K \B<£l x (t) - B™ (t) | < D m + 2M m+1 , (31) 

except for an event of probability at most (8 + 2 3 ~ 2C )(if 2 2m ) 1 ~ c . 

Thus we have to find an upper estimate M m . For that the large deviation 
inequality (fT2"]) will be used. By 0, the increment of Bm(t) on [t kl t k+ \] is 



A m , k = \B£\t k+1 )-B£\t k )\ 



rfg+ n E [(k + l-r) H ^-(k-r) H -^]x m (r+l). 

V 2/ r— — oo 
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Then a similar argument can be used as in the proof of Lemma 21 see e.g. part 
(a) there: 



T 1^3-W S r{A m , k ) = £ [(fc + i- P )*-*-(*- P )*-i 



r— — oo 
fc-2 



l+(2^-*-l) + £ (k 



T— — O0 

k-2 



< 1+ (2^-4 - l) + {k 



\2H-1 



\2H-1 



k — r 



1 



H - - 



H— 7 



1 



< 1+ H 



H 



1 



2 J (k - r) 2 

2 



2 2-2H ~ 2\ 2 



Hence taking N = K2 2m and C > 1 in (fT2"|) . and using (TH))) too, one obtains for 
m > 1 that 



M m = max |An,fc| 

Kfc<if2 2m 



5/V2 



r(F + i) 



2 



(1 -fO^CClog.-fiO* mi2- 2Hm , (32) 



with the exception of a set of probability at most 2 (i"C2 2m ) , where X > 
and C > 1 are arbitrary. 

Then substituting this and Lemma 0] into (f3~T|) . it follows that when K > 0, 
C > 3, and m > 7714(C), 

- < a{H,K)m2-^ m (33) 

except for an event of probability at most 8.125(X2 2m ) 1_c ' where a(H,K) and 
P(H) are the same as in Lemma E] Remember that the rate of convergence in 
(|32"1) . just like in parts (a) and (c) of the proof of Lemma |4j is faster than the 
one in parts (b) and (d) of that proof. Apart from constant multipliers, the 
result of (f3"2"j) has the same form as the results of (a) and (c) there. Since in the 
statement of this theorem we simply replaced the faster converging factors by 
the slower converging ones, the constant multipliers of (13"2l can be ignored if m 
is large enough. This is why the a(H,K) defined by Lemma H] is suitable here 
too. 

In the second part of the proof we compare Bm (t) to B m _^j(t), where j > 1 
is an arbitrary integer. If K > 0, C > 3, and m > 1714 (C), then (13"3"|) implies 
that 

m-\-j— 1 

m^ K \B^ 3 (t)-B^(t)\< J2 Q ^ K \B[ H + \(t)-B ( "\t)\ 



< Y a(H,K)k2-^ k < a(H ' K) 9 m2-*-""" 
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Hence one can get that 

P {sup max \B^Ut) B^(t)\ > / a ^ K) 2 m2~^A 

OO 

< 8 - 125 (^2 2fe ) 1_C < 9 (K2 2m )^ C . (34) 

By the Borel-Cantelli lemma this implies that with probability 1, the sample 
paths of B ( " ) (t) converge uniformly to a process W^ H ^ (t) on any compact inter- 
val [0,-K - ]. Then W^ H '(t) has continuous sample paths, and inherits the proper- 
ties of Bm (t) described in Section 3: it is a centered, self-similar process with 
stationary increments. As Lemma[S]below implies, the process {W^ H > (t) : t > 0) 
so defined is Gaussian. Therefore (t) is a fBM and by (|3"4")l the convergence 
rate of the approximation is the one stated in the theorem. □ □ 

The aim of the next lemma to show that integration by parts is essentially 
valid for @ representing W^ H '(t), resulting a formula similar to (fTUf. Then 
it follows that (W^ (t):t> 0) can be stochastically arbitrarily well approxi- 
mated by a linear transform of the Gaussian process (W(t) : t > 0), so it is also 
Gaussian. 

Lemma 5. LetWW(t) be the process whose existence is proved in Theorem [H 
above for H 6 (i, 1), or, by a modified construction, in Theorem^ below for 
any H G (0, 1). Then for any t > and e > there exists a So > such that 
for any < 5 < Sq we have 



> e < e, (35) 



where 

wf\t):=-[ h' s (s,t)W(s)ds, (36) 

J[-$,-5\U[0,t-S\ 

and h(s,t) is defined by 0). (W$ (t) is almost surely well-defined pathwise as 
an integral of a continuous function.) 

The lemma shows that as S — > 0+, W$ (t) stochastically converges to 
W( H \t) when H>\, while W^ H) (t) has a singularity given by the extra term 
in \35)) when H < |. (If H = | then Wg H \t) = and the lemma becomes 
trivial.) 

Proof. Fix t > and e > and take any S, < 8 < t. Let us introduce the 
notation, cf. ([8]): 

B%{t lm) )= £ Mtr ' AM( ^'^^ M) ^(^)At, (37) 

t r ei m ,s 

where I m j = \ (^) (m) < U < S( m ) U (0 < t r < t (m) - <5 (m) ] and the abbre- 

S(m) 



'(m) 

viation S( m \ = [s2 2m \2~ 2m is used for s = t,S, and —1/5 (an empty sum being 
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zero by convention. Then we get the inequality 



W {H) {t)- wl H \t) 



(H)l . 5 H -?W(t-5) 



<\W^ H) (t)-Bif\t {m) )\ 



n( H )(t \ R (H ht ) \ m) 2 B m (t (m) - 8 {m) + At) 



+\B%i(t {m) )-wr>(t)\ 



+ 



6 (m) 5 B ™ (*(m) - 6 (m) + At) gH- h\\r(t - 5) 



Y{H+\) 



Y{H + \) 



(38) 



First we have to estimate the second term on the right hand side as 5 — > 0+, 
uniformly in m (this requires the longest computation): 



where 



E 

*(m) <*r<t(m) 



T{H+\) 

h(t r - At,f( m )) - h(t r ,t( m )) 
At 



B m (t r )At 



6 (m) 2 B m{t(m) ~ <*(m) + At) 



r(F + J) 



x - /i(t r - At,t (m) ) -h{t r ,t (m) ) 

fm,6= 2^ a7 D m \t r ) l\t 



and 



— <5(m) <*r<0 

-oo<t r <(^-) (m) 



At 



/t(t r - At,t( m )) - /t(t r ,t( m )) 
At 



B m (t r )Ai. 



Then "summation by parts" shows that 

E m ,s = ^2 h(t r ,t( m -))[B m (t r+ i) — B m (t r )]. 

t(m)-*(m)<*r<t(m) 

(This is the point where the extra term in the definition of £? TOj( 5 is needed.) 
Thus 



var (r (fr + 1 ] /:,....- 



= t 



2ff-l 

(m) 



E 



t(m)-«(m)<tr<t(m) 

2H-1 



2m 



t(m)-«(m)<*r<t(m) 



1 - 



t(m) 



At 



< t 



2H- 
(m) 



<tf" 5 2H 
("0 < ° 

2H ~ 2H : 
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for any to > 0. Then by the large deviation inequality (jlip . for any m > and 
for any C > 0, 



P<^ |# m ,i| > (2Clog, 



r(/r + §)(2//)4 



< 25 c 



(39) 



Similarly as above, the definition of F m g can be rewritten using "summation 
by parts" that gives 

F m ,s = £ h(t r ,t( m ))[B m {t r+ x)— B m (t r )]+h(— (y( m ),t( TO ))-B m (— 6( m )+At). 

-6( m) <t r <0 

The definition of F m ^$ shows that it is equal to zero whenever 8 < At, therefore 
when giving an upper bound for its variance it can be assumed that S > At. 
Thus 



Var i I i 11 -\) F„,, 



£ [(t (m) +t v ) H -i-£ h 



0<t„<(5 ( „ 

x {5 {m) - At) 

0<t„<5 (m) 



At- 



(t(m)+S im) ) H k -t" m) h 



t{m) 



2H-1 



f v 



IH-l 



At 



(t(m) +$( m )) 2H 1 + S f m ) 1 



>(m) 



< t 



2H-1 

(rn) 

+ 2H 
2H 



Hra) 



2H-1 



-(m) 



2H-1 



dt + 2t lni) 1(5 (m) + 5 (m) 



2H 



r \2H / c x 2ff 



c (m) 



5(m) 



, ,2ff-lr , r2H 

+ 2 *(m) 



/ 6 ,2H-lc , (m) , ,2H-1 r , c2/f ^ ',2H-lc, ° c2H 
^ 2 *(m) <V) + 2ff + 2t (m) d (™) + V) < g* d + 



7 



_3_ 

2H~ 



So by the large deviation inequality (jlip . for any to > and for any C > 0, 



|F m , s | > (2Clog, (- 



r(i? + i) 



< 26 L 



(40) 



Proceeding in a similar way with G m> s, one obtains that 



G m! a — £^ /i(t r , t( m ))[S TO (t r _|_i) — B m (t r )] 

_ 00<tr< (^i) (m) 
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Hence 



< 



Var(r (H+-)G„ hS )= Y 



(t{m) + tv) H 2 — t v 



At 



(T i ) (m) <t " <00 



(*(™) ( S )(m) 

E C- 1 

-(-r) (m) <*««» 

/ \ 2H-1 

+ (~ (~r)(m) 



-(-&)t 
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At 



H "2 



C(m) 
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<5 /(m) 
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H 



< H 



<-(m) 



E e- 3 A*+(-(^) (m) 



2H-2 



,2H-3 



(-(^) 



(m) 



2H-2 



3(H 2^ ,2 A - J / / 

2(1- V 



So again by the large deviation inequality (|lip . for any m > and for any 
C > 0, 



P< |G m , 4 | > 2Clog 



1\ \ = i? 



rffr- 



i -# 



ts 



l-H 



< 28°. (41) 



Combining (|40p and (|4ip . it follows that there exists a 8q > such that 
for any < 8 < 8q and for any m > 0, 



flW (t(ro) )_ B w (t(fn)) 



if— - 

<^(m) 2 ^m(*(m) _ <^(m) + At) 



> - > < 



After the second term on the right hand side of (|38|) we turn to the third 
term. Take now any 8 £ (0, 8q) Since h(s,t) has continuous partial derivative 
w.r.t. s on the intervals [—1/8, —8] and [8,t — 8] and by Theorem [TJ B m a.s. 
uniformly converges to the Wiener process W on these intervals, comparing (|5o]l 
and ([37JI shows that with this 8 there exists an m such that 



{\B { "l{t [m] )-Wf\t)\>-^ 



< 



Theorem [T] also implies that m can be chosen so that for the fourth term in 
one similarly has 



5 {m) 2 S m(*(m) _ 8 (m) + Ai ) <5 ff ^VF(t - 8) 



T(H+\) 



> - > < 
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Finally, Theorem [5] (or, with a modified construction, Theorem [3] below) guar- 
antees that m can be chosen so that the first term in (|3"5|) satisfies the same 
inequality: 

p{|WW(t)- B^{t)\ >{}<- 4 - 
The last four formulae together prove the lemma. □ □ 

5 Improved construction using the KMT approx- 
imation 

Parts (b) and (d) of the proof of Lemma [U gave worse rate of convergence than 
parts (a) and (c), in which the rates can be conjectured to be best possible. 
The reason for this is clearly the relatively weaker convergence rate of the RW 
approximation of ordinary BM, that was used in parts (b) and (d), but not in 
parts (a) and (c). It is also clear from there that using the best possible KMT 
approximation instead would eliminate this weakness and would give hopefully 
the best possible rate here too. The price one has to pay for this is the intri- 
cate and "future-dependent" procedure by which the KMT method constructs 
suitable approximating RWs from BM. 

The result we need from Komlos, Major, and Tusnady (1975, 1976) is as 
follows. Suppose that one wants to define an i.i.d. sequence X\, X%, . . . of 
random variables with a given distribution so that the partial sums are as close 
to BM as possible. Assume that E(Xfe) = 0, Var(X k ) = 1 and the moment 
generating function E(e" Xfc ) < oo for \u\ < uq, uo > 0. Let S(k) — X\ + 
■■■ + X k , k > 1 be the partial sums. If BM W(t) (t > 0) is given, then 
for any n > 1 there exists a sequence of conditional quantile transformations 
applied to 1F(1), W(2), . . . , W{n) so that one obtains the desired partial sums 
5(1), 5(2), . . . , S{n) and the difference between the two sequences is the smallest 
possible: 

pj max \S{k)-W(k)\> C \ogn + x\ < K e- Xx , (42) 

[0<k<n J 

for any x > 0, where Cq,Kq,X are positive constants that may depend on the 
distribution of X kl but not on n or i. Moreover, A can be made arbitrarily large 
by choosing a large enough Co- Taking x = Cq logn here one obtains 

p{ max \S(k)-W(k)\ >2C lognl < K Q n- XCo , (43) 

I 0<k<n J 

where n > 1 is arbitrary. 

Fix an integer m > 0, and introduce the same notations as in previous 
sections: At = 2~ 2m , t x = xAt. Then multiply the inner inequality in (|4"3")l 
by 2~ m and use self-similarity ((T|) of BM (with H = I) to obtain a shrunken 
RW B* m (t k ) = 2- m S m (k) (0 < k < K2 2m ) from the corresponding dyadic 
values W(tk) (0 < k < K2 2m ) of BM by a sequence of conditional quantile 
transformations so that 

max \B* m {t k )-W{t k )\ < 2C^T m log, (K2 2m ) < 5C log, K m2" m , (44) 
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with the exception of a set of probability smaller than Ko(K2 2m )~ XCo , for any 
m > 1 and K > 0. [Here (fT9"]l was used too.] Then (|4"4l implies for the difference 
of two consecutive approximations that 



P max \B* m+1 (t k ) - B* m (t k )\ > 10C log, K m2~ 

< 2K Q (K2 2m )~ XCa (45) 

for any m > 1 and K > 0. This is exactly that we need to improve the rates of 
convergence in parts (b) and (d) of Lemma SJ 

Substitute these KMT approximations B* n (t r ) into definition ([7]) or ([5]) of 
Bm (tk ) ■ This way one can obtain faster converging approximations Bm^ of 
fBM. Then everything above in Sections 3 and 4 are still valid, except that one 
can use the improved formula P5"|) instead of Lemma [3] at parts (b) and (d) in 
the proof of Lemma [U This way, instead of pi"]) one gets 

f 23C log, K m2- 2Hm ifO<H<±, 
i<S- lYmM ~ { 15C log t KK H -im2-™ if \ < H < 1, (4bj 

for any m > 1, except for a set of probability smaller than 2Ko(K2 2m )~ xc ° . 
Also by (H5]), instead of (|2"5|) and ([2"B]) one has the improved inequalities: 



max 

{j-\)L<U<jL 



\B* m+1 {-t v ) - B* m (—t v )\ 



f 10C log,Lm2-" 1 ifj = l, 

" \ 14C log t j log,Lm2- m if j > 2, 

with the exception of a set of probability smaller than 2Ko(jL2 2m )~ xc ° , where 
m > 1. If Co is chosen large enough so that ACo > 2, then (|4"?| holds simulta- 
neously for all j = 1, 2, 3, . . . except for a set of probability smaller than 

oo 

2K (L2 2m )- XCo J2r XCo < -^K (K2 2m y XCo . (48) 
i=i 

(Remember that we chose L = AK in part (d) of the proof of Lemma |U) Then 
using this in part (d) of Lemma 21 instead of (|2"T|) one needs the estimate 



oo 
J=2 



/OO 
x H ~i log,a;dx < 1 (0 < H < 1) 



Then instead of (|28l) and (j2"9")l . the improved results are as follows. First, in the 



case ^ < H < 1 one has 



max \U m k\ 

Kk<K2 2m 



< 10C \og,Lm2- m h H -i + 14C log,Lm2-"^ 



H- 1 - 
2 



< 









^18 + 112 




) 









C \og*KK H -tm2~ m (49) 
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for any m > 1 and Co large enough so that XCq > 2, except for a set of 
probability smaller than given by (|4"5)l . 

Now in the case < H < 1 it follows that 



max 

Kk<K2 2 " 



\u, 



m,k | 



< WC \og,Lm2~ m -(At) H -i 



UC log*Lm2- m - 



H -2 



H- 



< C log, A to 36-2 



-2Hm 



H 



(50) 



for any to > 1 and Co large enough so that ACo > 2, except for a set of 
probability smaller than given by (|4"5)l . 

As a result, there is convergence for any H £ (0, 1). Since the KMT approx- 
imation itself has best possible rate for approximating ordinary BM by RW, it 
can be conjectured that the resulting convergence rates in the next lemma and 
theorem are also best possible (apart from constant multipliers) for approximat- 
ing fBM by moving averages of a RW. 

Lemma 6. For any H £ (0, 1), m > 1, K > 0, C > 1, and Co large enough, 
we have 



where t k = k2- 2m 

(log* 3 



max \B^l(t k ) - B< H Ht k )\ > «*m2^'W™ 

0<t k <K + 

< 4(A2 2m ) 1 - c + 3A (A2 2m )- ACa , 
/3*(H) = min(2i?,l), a* = a*(H, K,C,C ), 



a ' 



if He (0,i) , 



10C5 



(l-H)i 



+ Co (log, K) 5 59 + 84 



A' 



ff- 



r h 



10C*- 



Iff 



(1-^4 



CoOog.Jf) 3 33 + 112 



H 



A 



i/i? G (5, l), and i/ie constants X, Co and Ao are defined by the KMT approx- 
imation ^Sfy with Co is chosen so large that XCq > 2. [The case H — | is 
described by 



Proof. Combine the results of parts (a) and (c) in the proof of Lemma 0] and 
the improved inequalities above, that is, apply (fI5]l. ([2"0]). (gBJ), (JUJ), (gSJl, and 
(|50p . Here too, we simply replace the faster converging factors by the slower 
converging ones, but the constant multipliers of faster converging terms cannot 
be ignored, since the lemma is stated for any m > 1. □ □ 

Now we can extend the improved approximations of fBM to real arguments 
by linear interpolation, in the same way as we did with the original approxima- 
tions, see pop. This way we get continuous parameter approximations Bm (t), 
(t > 0) for m = 0, 1, 2 } . . . , with continuous, piecewise linear sample paths. Now 
we can state the second main result of this paper. 
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Theorem 3. For any H € (0, 1), the sequence B^ H) (t) (t > 0, m = 0, 1, 2, . . . ) 
a.s. uniformly converges to a fBM W^ H \t) (t > 0) on any compact interval 
[0,K], K > 0. If m > 1, K > 0, C > 2, and Co is large enough, it follows that 

p { og?K |w«(.) - itfO(„| > _|_ im2 -»- ( », J 

2m\l-C , ^ / 7^n2m\-ACo 



< eo^™) 1 +4if (^2" m ) 



with 



* 10 ^.i /i 

" + F(flTl) Cl(log - A) 



-4 



where a* and f3*(H) are the same as in Lemma\^ (In other words, in the 
definition of a* in Lemma\fi\the constant multiplier 10 has to be changed to 20 
here.) The constants A, Co, Kq are defined by the KMT approximation \4'2fy with 
Co is chosen so large that ACo > 2. [The case H — \ is described by j44\ )-l 

Proof. The proof can follow the line of the proof of Theorem [2] with one ex- 
ception: the constant multipliers in (|32[) and consequently in pip cannot be 
ignored here. This is why the multiplier a* of Lemma [HI had to be modified in 
the statement of the theorem. □ □ 

It can be conjectured that the best rate of approximation of fBM by moving 
averages of simple RW's is 0(N~ H logiV), where N is the number of points 
considered. Though it seems quite possible that definition of B*n H%> {t) above, 
see ([7]) with the KMT approximations B^(t r ), supplies this rate of convergence 
for any H £ (0,1), but in Theorem [3] we were able to prove this rate only 
when H g (0, \). A possible explanation could be that in parts (b) and (d) of 
Lemma [4] we separated the maxima of the kernel and the "integrator" parts. 

As a result, the convergence rate we were able to prove when | < H < 1 
is the same 0(N~^ log TV) that the original KMT approximation (|44|) gives for 
ordinary BM, where N = K2 2m , though i n this case the sample paths of fBM 



are smoother than that of BM. (See e.g. Decreusefond and Ustiinel (1998) .) 
On the other hand, the obtained convergence rate is worse than this, but still 
thought to be the best possible, 0(N~ H log N), when < H < i, which heuris- 
tically can be explained by the more zigzagged sample paths of fBM in this 
case. 
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